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Abstract 

We study theoretically the viscoelastic properties of sheared binary fluids that have strong dynamical 
asymmetry between the two components. The dynamical asymmetry arises due to asymmetry between 
the viscoelastic stresses, particularly the bulk stress. Our calculations are based on the two-fluid model 
that incorporates the asymmetric stress distribution. We simulate the phase separation process under an 
externally imposed shear and compare the asymmetric case with the usual phase separation under a shear 
flow without viscoelastic effects. We also simulate the behavior of phase separated stable morphologies 
under applied shear and compute the stress relaxation. 

PACS numbers: 83.60,83.50.Ax,64.75.+9 
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I. INTRODUCTION 



Many fluids can exhibit viscoelastic behavior, i.e. their response to deformation is intermediate 
between that of solids and fluids. For short times, the response is elastic and the stress is propor- 
tional to the applied strain. On the other hand, in the long time limit, a fluid like response with 
stress proportional to the strain rate is observed. The effect of viscoelasticity on the morphologies 
of phase separating polymer solutions and blends has been recently experimentally studied by 
TanakaM. It is now well established that due to an asymmetry between the viscoelastic properties 
of the two phases, a transient network of the more viscoelastic phase is observed. While there 
is some understanding of the influence of viscoelastic effects on phase separation, little is known 
on how viscoelasticity effects phase separation when an external shear flow is applied. Although 
effects of shear on newtonian (purely hydrodynamic) binary fluids have been well studied^ "^ and 
dynamical steady states with stringy morphologies have been observed, it is also important to 
study similar situations for binary fluids where viscoelastic effects can dominate. 

Theoretically, the coupling between the stress and concentration fluctuations was investigated 
by Helfand and Fredrickson^ in the context of shear flow in polymer solutions. Doi and Onuki^ 
studied this coupling by introducing a two-fluid model that considers two different velocities for 
the two components. Taniguchi and Onukii have used the two-fluid model to study viscoelas- 
tic phase separation. However, a two-fluid model incorporating dynamical asymmetry due to 
asymmetric distribution of bulk stresse s^-^i'"i" has successfully explained many features of the 
viscoelastic phase separation experiments of Tanaka, such as the transient network formation and 
eventual phase inversion. Recently, this theory has also been applied to study the viscoelastic 
phase separation in diblock copolymers'^. The two-fluid model was also used to simulate poly- 
meric foams under sheapi^. In this paper, we use this model to study phase separation under a 
shear flow for a binary fluid that is characterized by a strong asymmetry between the viscoelastic 
moduli of the two phases, particularly the bulk modulus. 

In addition to phase separation under shear flow, the behavior of stable phase separated states 
under a deformation is also important. Recently, Buxton and Balazas^ have studied the mechan- 
ical properties of phase separated morphologies using a block and spring model. This model can 
only model elastic properties and cannot describe the viscoelastic response. Thus we study the be- 
havior of fully phase separated stable states under a constant shear based on the two-fluid model. 
The time evolution of mechanical stresses generated for such phase separated states is also studied. 
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The paper is organized as follows. In Sec. 2, we discuss the two-fluid model used in our 
simulations. Section 3 describes results on viscoelastic phase separation under shear flow while 
section 4 is devoted to the response of stable phase separated morphologies subjected to steady 
shear. We conclude the paper with a summary and discussion in Sec. 4. 



n. TWO-FLUID MODEL 



The coupling between viscoelastic stresses and diffusion is studied using the two-fluid model. 
We consider a dynamically asymmetric mixture of a viscoelastic fluid (A) and a purely hydrody- 
namic (non- viscoelastic) fluid (B). In this model, different velocity fields, VAir,t) for the vis- 
coelastic phase and t'B(r, t) for the non-viscoelastic fluid are introduced. The average velocity is 
given by Vaverage{r, t) = (f){f, t)vA{r., t) + (1 — 0(r, t))vi{f, t), where 0(r, t) is the concentration 
of the viscoelastic fluid. The free energy of mixing can be written in terms of the concentration 
0(r, t) as 

Frm. = j dr[ct>\l - 0)' + (C/2)(V0)2], (1) 

where C is the concentration gradient coefficient that is assumed to be constant in the present 
work. The free energy Fmix is a simplified version of the Flory-Huggins^^ type free energy. The 
applied shear flow is implemented by considering an average velocity Vaverage = VappUed + v, 
where v is the contribution to the average velocity due to concentration fluctuations and Vappiied is 
the contribution due to the external shear flow given as, 

V applied = s\x — XqU . (2) 

This represents an external shear (symmetric about Xq) applied along the y direction and s is the 
shear rate. The appropriate equations of motion for incompressible viscoelastic binary fluids under 
external shear flow are 

86 - , ,dd) ^ 0(1-0)2 ^ ^ ^ 

^ + ^. V0 + s X - xohr = V ■ . V -U-V-a]. (3) 
dt dy C 

<-> <-> ^ 

Here cr represents the viscoelastic stress tensor. The osmotic tensor n is defined as V- 11= 

0V{20(1 — (f)Y — 2(1 — 0)0^ — CV^0}. The quantity ( is the friction coefficient associated 

with the velocity tr^ — f^. The dynamics of the average velocity field v is given by the usual 

Navier- Stokes equation 

dv ^ ^ ^ ^ ^ „ 
p-^ = -V- n +VP + V- a +r/,VV, (4) 



where a pressure has been introduced to account for the incompressibility constraint. If we use the 
overdamped limit || = 0, then u 
be expressed in Fourier space as 



overdamped limit || = 0, then under the incompressibility condition V ■ ?7 = 0, the velocity can 



Vk = Tk ■ (-V- n +V- 

where {)k represents the Fourier transform and Tk is a tensor defined in fourier space as 

1 



(5) 



Tk 



kk 
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where J represents the unit tensor. Next, we have to specify a constitutive law for the viscoelastic 
stresses. We use the Maxwell model according to which the equation of motion for the stresses is 
given by 



dt 
d a. 



M, 



ms{4>) M, 



(7) 



dt r,(0) 

where and are the internal molecular relaxation times associated with the bulk and shear 
stresses respectively. The matrix M is given as 



M-- 



dx 

^ dy dx ' 



^ dy dx > 

C2 p-VAy 

dy 



(8) 



The velocity v~X represents the velocity of the viscoelastic phase defined as 



V^A = V + Vapplied " 



'1- 



V- n -V- a 



(9) 



The above definition of M assumes that the stress acts on the viscoelastic fluid only. The final 
stresses (T;, and cr^ are calculated as ct;,' = {\/2)Tr{ai)) I and cr<, =(7^ —{\/2)Tr{as) /. The 
total stress is calculated as a=ai, + . The model decribed in equations (1) to (9) can be used 
to simulate the viscoelastic phase separation under shear flow with shear rate s. 



Ill, VISCOELASTIC PHASE SEPARATION UNDER SHEAR FLOW 



We describe our results on simulations of viscoelastic phase separation, both with and without 
shear. The model described in Sec. II is discretized on a 128 x 128 grid with periodic boundary 
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conditions, using a spectral approach for the velocity equation and explicit, central difference 
schemes (Ax = 1, At = 0.005) for the rest. 

Following Tanaka and Araki^iS, we choose rrih = mh^9((f) — 0o) (0o is the initial concentration 
before quenching), = mj^ct)^, r^, = Th^cf^ and = Ts^cf^. With this choice of the moduli 
and time scales, Tanaka and Araki were able to reproduce a number of experimental features of 
viscoelastic phase separation^. 

We first simulate phase separation without an external shear(s = 0). We consider a 50 — 50 mix 
with initial concentration fixed at 0o = 0.5. The field 0(r, 0) is initialized by small fluctuations 
of order 0.001 around 0o and the model in Sec. II is simulated for three different cases. Case I 
corresponds to hydrodynamic phase separation in the absence of viscoelastic effects. Here we set 
all the stresses to zero and solve Eqs. (3) and (5) with only one velocity field v. Case II is the full 
asymmetric viscoelastic model with rfs = \, C, = 0.1, m}P = 5, m^" = 0.5, Tf," = 10 and = 50. 
To clarify the role of the bulk stress, in case III we switch off the bulk modulus, i.e., m^^ = 0. 

Figure 1 shows the phase separation process for the three cases. The left column shows the 
non-viscoelastic (hydrodynamic) case, middle column depicts the full viscoelastic model with 
asymmetric bulk modulus and the right column represents the phase separation for viscoelastic 
case without bulk modulus. This figure shows the importance of the bulk modulus in the formation 
of a network of the more viscoelastic component. For the case with the bulk modulus, initially, 
holes of the non viscoelastic phase are formed in a partially phase separated matrix of the more 
viscoelastic phase. As these holes grow, the area of the viscoelastic region decreases and a network 
like structure is formed. These results are consistent with earlier results of Tanaka and Araki^. For 
the full viscoelastic model, we do not observe the breaking of the network within the time interval 
we simulated. The case without the bulk modulus appears to be similar to the purely hydrodynamic 
case and there is no network formation. Figure 1 also shows that viscoelastic effects slow down 
the domain growth, as can be inferred by comparing the domain patterns at time t = 1000 for all 
the three cases. 

We have also monitored the evolution of the viscoelastic stresses during the phase separation. 
Figure 2 shows the behavior of the average normal stress o-xx—Uyy, the average bulk stress (Txx+o'yy 
and the average shear stress a^y for the case with bulk stress as well as that for the case without 
bulk stresses. The normal and shear stresses do not decay completely in the time simulated and 
exhibit small fluctuations of order 0.001, possibly due to the motion of the viscoelastic phase. The 
behavior of the bulk stress is worth noticing as for the case with a bulk modulus, the early stage 
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of phase separation is associated with significant compressive stresses. The compressive stresses 
are associated with the formation of the network of the more viscoelastic phase. The network of 
partially phase separated viscoelastic phase shrinks as the phase separation proceeds resulting in 
the compressive stresses observed in the early stages. 

In Fig. 3, we display the phase separation process when an external shear with shear rate 
s = 0.01 is applied. As in the case without shear, we consider three cases corresponding to the 
hydrodynamic case without viscoelastic effects, full viscoelastic model with bulk modulus and vis- 
coelastic model without bulk modulus. At this shear rate, the very early stages are not influenced 
significantly by the applied shear for all the cases. At later times, there is an underlying tendency 
for the domains to align along the shear direction, although the alignment is not complete within 
the time simulated {t = 2500). The case with the bulk modulus exhibits network morphologies in 
the earlier stages even for this case. The alignment is also slower in the initial stages for the case 
with the bulk modulus. 

The time evolution of viscoelastic stresses for this value of shear is shown in figure 4 for the 
case with bulk stress as well as case without bulk stress. Fluctuations due to random motion of 
domains in the early stages are observed for both normal and shear stresses and the amplitude of 
the fluctuations decreases at late times as the domains tend to align along the flow direction. For 
this case also, the bulk stress for the case with bulk modulus becomes negative in the early stages, 
indicating the shrinking of the network of the viscoelastic phase. 

Finally, we consider phase separation under a relatively large shear rate s = 0.1. Figure 5 
shows the phase separation at this shear rate for all three cases. For this case, the external shear 
effects dominate the phase separation and the tendency for the domains to align along the shear 
to form stringy patterns is much stronger than the earlier case. The purely hydrodynamic model 
rapidly forms the stringy phase, whereas for the viscoelastic cases, the intermediate stages are 
characterized by complex motion of the domains. For the case with the bulk modulus, the early 
stages do not show the formation of a well developed network as the shear dominates over the 
internal dynamics. The snapshots at t = 2500 show that a stringy phase seems to form at long 
times for the viscoelastic cases also, although the width of the lamellae is larger compared to the 
hydrodynamic case. 

Figure 6 shows the time evolution of the viscoelastic stresses for this case. After transient 
fluctuations, the normal stresses decay to zero as the steady state is established. Compressive bulk 
stresses are observed even for this case but the magnitude is smaller for this case, compared to 
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the earlier cases. This is due to the fact that the shear flow is so fast that the transient network 
of the viscoelastic phase is not allowed to fully develop. This is clear from the snapshot at time 
t = 50 for the viscoelastic case. The behavior of the shear stresses is interesting. At long times, 
the shear stresses saturate to non zero values. The saturation stresses are related to the effective 
shear viscosities for this binary fluid. 

IV. PHASE SEPARATED STABLE STATES UNDER SHEAR 

So far we have considered the effects of shear on binary fluids undergoing phase separation. 
In this section we apply shear on stable phase separated morphologies. The aim of this section 
is to explore the possibilities of using this model to study mechanical properties of polymer solu- 
tions and blends. Recently, Buxton and Balazs'"^ have used a block and spring model to simulate 
deformation of a random two-phase morphology obtained from a Cahn-Hilliard simulation. In 
this approach, there is no coupling between the morphologies and the deformation. Such an ap- 
proach can only be used to study linear elastic behavior and does not describe the full viscoelastic 
response. The morphologies cannot evolve in response to deformation. However, the two-fluid 
framework used in this paper explicitly incorporates a coupling between the stress and the concen- 
tration and can describe deformation induced morphological evolution. 

Here, we consider a simple phase separated morphology below the coexistence temperature. 
We consider a 50 — 50 (0o = 0.5) binary fluid that exists as a band of the fluid A sandwiched be- 
tween layers of fluid B. This is a stable phase separated configuration. Polymer blends can show 
such macro-phase separation in equilibrium and block copolymers also micro-phase separate into 
lamellar morphologies^^. We apply shear in the same way as discussed in the previous section. We 
consider two different shearing conditions. In one case the shear is normal to the A/ B interface 
and in the other case, the shear is parallel to the interface. Such a situation (for block copolymers) 
has been investigated experimentally^ and theoretically using hydrodynamic models^^. The mor- 
phological evolution for the case of shear normal to the A/ B interface is depicted in Fig. 7. All 
three cases that have been discussed in the previous section are simulated for this case also. For all 
the cases, initially the interface starts to move in response to the shear and at long times there is a 
tendency to form many bands aligned along the shearing direction. The bands appear to be wider 
for the viscoelastic cases in comparison to the hydrodynamic case. However, we do not find much 
difference between the case with bulk modulus and the case without bulk modulus. Interestingly, 
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for the case when the shear is applied parallel to the interface, the band remains stable and no mor- 
phological evolution is observed for all the three cases. This is consistent with the shear induced 
reorientation of lamella for block copolymers observed in experimentSiii and theoretical models^^. 

Figure 8 shows the evolution of stresses for the full viscoelastic model (with bulk modulus) for 
both cases (shear parallel and perpendicular to the interface). For shear parallel to the interface, 
the normal and bulk stresses remain zero during the shearing process. This is due to the fact that 
there are no morphological changes and only shear stresses are generated. For the case of shear 
normal to the interface, normal as well as bulk stresses are generated due to shear induced motion 
of the interface depicted in Fig. 7. For a short interval, weak compressive stresses are observed. 
The early time behavior of the shear stress also differs from the case when the shear is applied 
parallel to the interface. This figure demonstrates the role played by the morphological evolution 
on the viscoelastic properties of binary fluids. For the case without bulk modulus, very similar 
stress evolution (not shown here) is observed. However, no bulk stress is generated. 

Finally, we study the relaxation of the morphologies after the applied shear is suddenly re- 
moved. We remove the shear at t = 50 corresponding to the morphologies shown in the bottom 
row of Fig. 7. Fig. 9 shows the final morphologies at t = 2150 after removing the shear. It is 
clear that the morphological relaxation is much faster for the non-viscoelastic case and the band 
reappears, although it is shifted from the original position. For the viscoelastic cases, the morpho- 
logical relaxation is very slow and the pattern does not evolve to the original band within the time 
span of the present simulations. 

V. SUMMARY 

We have investigated the viscoelastic properties of dynamically asymmetric binary fluids under 
shear flow. The dynamical asymmetry arises as the binary fluid constitutes one viscoelastic fluid 
and the other non-viscoelastic purely hydrodynamic fluid. We have simulated phase separation 
for this binary fluid both with and without shear. Only for the case with a nonzero bulk modulus, 
there is a tendency to form a network of the more viscoelastic phase in the initial stages of phase 
separation, both with and without shear flow. However, for the high shear rate, the network phase 
is short lived as the shear has a tendency to suppress network formation. For the high shear rate, 
in the long time limit, a stringy phase is observed even for the viscoelastic cases, however the 
length scales associated with the stringy phase are larger compared to the purely hydrodynamic 
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case. Thus, viscoelastic effects enhance the extent of phase separation under shear. We have also 
studied the temporal evolution of the effective viscoelastic stresses during the phase separation 
process. Interfacial motion results in fluctuations of the shear and normal stresses in the early 
stages. Transient compressive stresses are also observed corresponding to shrinking of viscoelastic 
phase in the early stages. At long times, the shear stresses saturate to a nonzero value. This value 
is related to the effective viscosity of the two-phase fluid. 

We have investigated the effects of shear on stable phase separated morphologies. A lamellar 
structure of the binary fluid below the coexistence temperature is sheared, both normal and parallel 
to the interface. There are crucial differences between these two cases. For the shear parallel 
to the interface, no morphological evolution is observed while for the case of shear normal to 
the interface, the interface moves and splits into bands that tend to align along the shear. The 
evolution of effective stresses also depends on the direction of the interface. For example, for the 
shear parallel to the interface, no bulk and normal stresses are generated. On the other hand, for 
the shear normal to the interface, bulk, normal and shear stresses are generated due to interfacial 
motion. 

This paper demonstrates the role played by viscoelasticity on shearing of binary fluids. The 
simulations show that the underlying viscoelastic properties can significantly influence the shear- 
ing behavior. Further experiments are needed to test predictions of our simulations. The approach 
used in the present work can also be applied to other complex fluids such as block copolymers and 
microemulsions. 
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Figure Captions: 

1. Phase separation process without shear flow for the purely hydrodynamic case, viscoelastic 
case with bulk modulus and viscoelastic case without bulk modulus. 

2. Evolution of the normal, bulk and shear stresses for the situation depicted in Fig. 1. The 
stresses are calculated by averaging the local stresses over the entire system. 

3. Phase separation process with shear flow of shear rate s = 0.01 for the purely hydrodynamic 
case, viscoelastic case with bulk modulus and viscoelastic case without bulk modulus. Shear is 
applied along the y direction and periodic boundary conditions are assumed. 

4. Evolution of the normal, bulk and shear stresses for the situation depicted in Fig. 3. 

5. Phase separation process with shear flow of shear rate s = 0.1 for the purely hydrodynamic 
case, viscoelastic case with bulk modulus and viscoelastic case without bulk modulus. Shear is 
applied along the y direction and periodic boundary conditions are assumed. 

6. Evolution of the normal, bulk and shear stresses for the situation depicted in Fig. 5. 

7. Evolution of the band morphology under a shear flow of shear rate s = 0.1 for the purely hy- 
drodynamic case, viscoelastic case with bulk modulus and viscoelastic case without bulk modulus. 
Shear is applied along the y direction and periodic boundary conditions are assumed. 

8. Evolution of the normal, bulk and shear stresses for the case with nonzero bulk modulus 
(evolution depicted in the middle column of Fig. 7). 

9. Morphologies for the three cases at t = 2150 time steps after the shear flow has been 
removed. The shear flow was removed at t = 50 steps after the shear was applied (corresponding 
to the last row in Fig. 7). 
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